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RHIC is now confirming the discovery of the "New State of Matter", whose properties are grad- 
ually revealed by experimental data. It is not the free gas state, but well described by a fluid 
If^ . with very small viscosity. It is now highly desired to calculate the value of the viscosity from the 

\ fundamental theory, i.e., QCD just above the transition temperature. In this report we present our 

calculation of the transport coefficient of gluon system on 24^ x 8 lattice in the quench approxima- 
tion. Simulations are carried out in the range, lA <T /Tc < 24. In the temperature region slightly 
Q^' above the transition, where the perturbative calculation is not applicable, the shear viscosity(Tj) is 

^ • smaller than typical hadron masses. The bulk viscosity is consistent with zero within the range of 

error bars in 1.4 < T /Tc < 24. We compare our results with the perturbative calculations in large 
T /Tc region. It is found that the lattice and perturbative results are consistent with each other 
. there. The ratio rj /s is around 0.1 — 0.4 in T /Tc < 3 region and satisfies the KSS bound [p^t. In 

^ I order to estimate the contribution from high frequency part of the spectral function, we study the 

effects of a term p'"?^' proposed by Aarts and Resco[^. It is found that until the threshold mass 
becomes small, its effect is quite small, and that viscosity decreases as the threshold decreases. 
From these studies we think that although our result is obtained under an assumptions for the 
spectral function, it gives a reasonable estimation for ?](= ndp/d(0 at o = 0), and qualitative 
results will not be changed when the accurate spectral function is obtained. 
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1. Introduction 

The RHIC data are now revealing some properties of a "New State of Matter", wliicii is ex- 
pected to be quark-gluon plasma(QGP). Tiie jet quenciiing data suggest tiiat its mean free patii is 
short, and the phenomenological analysis of the elliptic flow shows that, it could be explained by a 
fluid model with small viscosity. Namely it is almost perfect fluid. As the viscosity is proportional 
to mean free path|jl|], this analysis also indicates that mean free path is short. The mean free path 
A is related to scattering cross section a and the number density « as A oc l/(a«). Therefore the 
small A, means that the "New State of Matter" is a strongly interacting system. 

It is a natural anticipation that transport coefficients of the QGP must be smaller than an or- 
dinary matter, like water or oil etc., because QCD coupling in the deconfined phase will be larger 
than the electromagnetic coupling constant 1/137, but its value must be calculated from funda- 
mental theory of QCD. Therefore it is now urgent to calculate transport coefficients by fully taking 
into account non-perturbative effects, because QGP near the transition temperature will be strongly 
interacting In Ref.[^ we published our results from lattice simulations in the 1.5 < T /Tc < 2.1 
region. In this report, the study is extended to higher temperature region, T /Tc < 24, and we com- 
pare our results with perturbative calculations there. We also discuss how our results depend on the 
spectral function of the green function in high frequency region[||]. 

On a lattice, calculation of the transport coefficients is formulated in the framework of the 
linear response theory [I,!]. 

T7 = - / dh' f dhe'^''-'^ r dt' < Ti2ix,t)Ti2ix'/) >ret (1-1) 



ri + i; = - [ d\' f dhe'^'^-'^ r dt' < Tn{x,t)Tn{x',t') (1.2) 



where 77 is shear viscosity, and l^, bulk viscosity. < T^yTpa >ret is a retarded Green's function of 
energy momentum tensor at a given temperature. In the quenched model, T^y is written by the field 
strength tensor as follows. 

r^v = 2Tr[F^aFva -- ^S^yFpaFpa] (1.3) 

The field strength tensor is defined by the plaquette operator, U^v{x) = exp {ia^gF^v{x))- 

The shear viscosity in Eq. |l.l| is also expressed by using a spectral function p of the retarded 
Green's function p(G))|j5|] as follows. 

r, = ;rlim^ = ;rhm^ (1.4) 

£0^0 (0 £0^0 d(0 

For the determination of the spectral function p(ft)), we use a well known fact that the spectral 
function of the retarded Green's function at temperature T is same as that of Matsubara-Green's 
function. Therefore our target is to calculate Matsubara-Green's function(G^ (f„)) on a lattice and 
determine p(co) from it. 
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2. Lattice calculation of 77 from Matsubara-Green's function Gp{t) 
Using the spectral function p{o)), Gp{t) is expressed as follows 

'^^ ' .l-ooi(0„-(0 Jo sinh{(0^/2) ' ' 

where (On = 2nn/p. Note that the last equality holds for the region < ? < P{Nt)- at f = and 
j3(= Nt) it could not be used. 

In order to determine the spectral function p{(o) from Gp{tn), which is given at discrete 
points ?„ = 1,2...,A^7' — 1 on a lattice, the most promising method may be the maximum entropy 
method(MEM). However to determine the p{(o) precisely by MEM, Gp{tn) on Nj > 30 lattice 
with high statistics is necessary. However, we realized that {tn) is very noisy [^], and we decided 
to start with a relatively smaller lattice of large statistics with some assumption for the spectral 
function p{(o). 

The simplest non trivial form is a Bright- Wigner type proposed by Karsch and Wyld|Q]. 

Ay Y 
pBW {(>>) = —{-, 0—-, ^ t)- (2-2) 

In this case, sheai" viscosity is given by T] = 4A 7^/(7^ + m^)^, similarly for bulk viscosity As 
this formula already has 3 parameters, in order to determine them by least squares, the lattice size 
in temperature direction(A^7-) must be A^r > 8. Then the minimum lattice size is 24^ x 8, to obtain 
a non trivial results. 

Simulations are carried out by using the Iwasaki's improved action and standard Wilson action. 
The simulation points are at j8 =3.05, 3.3, 4.5 and 5.5 for the improved action and j8=7.5 and 8.5 
for Wilson action. With roughly 10^ MC measurements at each j8, we determine the Matsubara- 
Green's functions Gj3(f„) as shown in Fig. [I| 

Although the errors of the Matsubara-Green's function are still not small in large t region, we 



can fit them with the spectral function p, given by Eq.p^ using a least square package SALS. The 
fitting range is 1 < f < 4. The viscosities are calculated from these parameters, and the errors are 
estimated by the jackknife method. The bin size of the jackknife analysis is taken to be 10^ and 
2 X 10^ for the improved and Wilson actions, respectively. 

The bulk viscosity is equal to zero within the range of error bars, while the shear viscosity 
remains finite. In order to determine the shear viscosity in physical unit, the lattice spacing <3(j8) 
must be determined, because, in the lattice calculation, rj x is obtained. The relation between 
lattice spacing a and p is studied phenomenologically in 2.2 < j8 < 3.8, and 5.58 < j3 < 6.5 re- 
gions for improved[|8|] and Wilson actions[^], respectively. In the larger j8 region, we adopt 2-loop 
renormalization group relation, 

«(« = X(^)"*«p(^), (23, 

with bo = 33/(487r2), and bi = (34/3)(3/167r2)2. The A is determined by the a(j3) at j3 =3.8 and 
6.5, and they result in T^/A = 1.26 and 47.35 for the improved and Wilson actions, respectively, 
where = 270MeV in the quenched model. The results for the shear viscosity are shown in Fig. 
^ with circle. It is seen that rj from both actions are consistent with each other. 
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Figure 1: Matsubara-Green's function Gnit) at j3 = 3.3 

3. Comparison with perturbative results 

As we have calculated rj at rather high temperature region, we compare our results with the 
perturbative results. The perturbative results are summarized as follows. The bulk viscosity ^ is 
zero ||5|, [I^. The shear viscosity in the next- to-leading log is given byljll]], 

where T71 = 27.126 and jl* /T = 2.765, and g is the running coupling constant. We use 2-loop 
formula for it, 

g-\T) = IhM^ij) + ^In (21n(M^)) (3-2) 
Note that the perturbative approximations are renormalization-scale /i dependent. The sim- 



plest choice would be = Tc/A.qcd ~ 1.03[|l2|]. However in order to see the ambiguity of the 



perturbative result due to /x, we change }X in the range 0.5 < /X < 2. The results are also shown 
in Fig. ^ with line. In the large T region, the /i dependence is not large. However at smallest T 
of each line, it shows a small increase. It means the break down of the perturbative calculation. 
However until the break down starts the agreement of the perturbative and lattice results are satis- 
factory. From this, we think that although our result depends on the assumption of the psw given 
in Eq. it may be a reasonable approximation for dp/dco at ft) = 0. 
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Figure 2: Shear viscosity in physical unit from lattice(circles) and perturbative calculations(hnes). The lines 
show the ambiguity due to the choice of the scale parameter in two-loop running coupling formula. 



Let us proceed to the study of the rj /s ratio, recently discussed in []l3|, The entropy on the 
lattice is reported by [||] and [ |l5| ] in T/r^ < 4.5, for improved and Wilson actions respectively. In 
the higher temperature region, where the data are not available, we use the high temperature limit 
values on A^r = 8 lattice. The results are shown in Fig. || with circles. For the perturbative result 
on entropy, we use the formula of hard thermal loop calculation up to order [ 12]. The results are 
also shown in Fig.^ with lines. 

As explained in the Fig.^, the slight increase at the small T region is observed for the per- 
turbative results, indicating breakdown of perturbative calculations. The agreement between these 
two calculations are not good enough. However we do not think that it is a real difficulty. Because 
if we decrease the scale parameter /i, the agreement is improved. An important information from 
these results is the qualitative magnitude of the ratio rj/s. In the T/Tc < 3 region, the ratio is 
small (0.1-0.4), where the perturbative results diverges. And for all the temperature regions we 
have studied, it is less than one. We don't think that it would become 10 times of present value, 
when the accurate determination of the spectral function is carried out. 

4. Discussions and conclusions 

4.1 Fit of Gp by other spectral functions 

Aarts and Resco have proposed an another form for the spectral function pp]. 
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Figure 3: tj /s obtained by lattice simulations(circles) and perturbative calculations(lines). The is a scale 
parameter in the running coupling constant of Eq.3.2 



p(co)=p""^'(«) + P^'^''(« 



bi+b2X^ + ... 

1 + C\X^ + C2X!^ + .. ' 



ft) 
T 



(4.1) 
(4.2) 



p%''(«) = e(ftJ-2m,,) 



J^(ftj2-4m2j5/2 



80712 w 



[«(ftj/2) + 0.5] 



(4.3) 



where cIa =N^ — 1, and n{(o) = 1/ (exp(ftj/r) — 1). We consider three parameters case, bi, ci and 
rrith. The fit by SALS could be made but c\ <0. 

We would like to remember that the spectral function p must satisfy the sum rule[p^], 

1 



-co 271 



(op{(o)d(o = 1. 



(4.4) 



In the full parameter representation, the cancellation of the divergence between p'""' and p'"^'' 
in high CO region results in the convergence of the integration. But in this three parameter form, 
the p could not satisfy the sum rule. Then some modification is necessary. But in this work we 



applied the unmodified form. Because the integration over (O in Eq. 2. 1 converges due to the term 
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cosh{(o{t — P /2))/sinh{(oP /2), except for f = and f = j8, we think that phenomenologically, a 
possible contribution from high frequency part of p(ft)) could be studied in this form. 

As the another trial, we study effects of p'"^'' for the shear viscosity rj. We assume p is given 



by p = Pbw + p^'^^, where psw is given by Eq. 2/2. By changing m,/,, the change of rj is studied at 



j8 = 3.3 for the improved action. When p'"^^'' is absent (m^/, = oo), 77a^=0.00225(201). If nith is set 
to be 5.0, 3.0 and 2.0, Tja^ becomes 0.00223(0.00191), 0.00194(0.00194) and 0.00126(0.00204), 
respectively. And at mti, = 1.8, the contribution from p^'^'^ becomes larger than the Gp{tn) of 
simulation at t = I, that fit could not be done. Generally as nith decreases the contribution from 
phigh increases and the p in the small (O region is suppressed. In this case, it results in the decrease 
of rj. 

In this case, the Pbw (without p'"^''), determined by SALS to fit the Matsubara Green's func- 
tion, satisfies the relation, 

■^(opBw{(o)d(0 < 1. (4.5) 

-c» 271 

As a full spectral function p will have contributions from scattering states etc., in addition to Bright- 
Wigner type, Eq j4.5| does not mean that the full spectral function violates the sum rule. 

4.2 Conclusion and further studies 

On 24^ X 8 lattice, we have calculated the Matsubara-Green's function and determine the trans- 
port coefficients of gluon plasma. The bulk viscosity is consistent with zero, while shear viscosity 
remains finite, within the error bars of present statistics. In the high temperature region, the agree- 
ment of the lattice and perturbative calculation is satisfactory for the shear viscosity. The lattice 
result on ratio rj /s in T /Tc < 3 is smaller than the extrapolation of the perturbative calculation and 
satisfies the KSS bound. 



Although our results depend on the form of the spectral function Pbw given by Eq.2.2, we think 
that the qualitative features will not be changed. Because as discussed in the previous subsection, 
our results are stable if the high frequency part of the spectral function p'"^'' is included. It is ex- 
pected that as p ( CO ) decreases faster than ft) and due to the factor , cosh {(o{t — ^ /2)) / sinh{(0^ /2) 
in Eq. |2.l| , the Matsubara Green's function Gp{t) will not be sensitive to the behavior of p(ftj) in 
high ft) region, except at f = and t = p. We think that the Tj and rj/s will not become 10 times of 
the present value when more accurate determination of the transport coefficients are carried out. 

However it is important to make a more accurate calculation of the transport coefficients, 
independent of the assumption of the spectral function. For that goal, we are starting the calculation 
of Gp{t„) on an anisotropic lattice, to apply the maximum entropy method. 
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